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The magnetic coupling between two concentric mesoscopic superconductors with non-zero thick- 
ness is studied using the nonlinear Ginzburg-Landau theory. We calculated the free energy, the 
expelled field, the total field profile, the Cooper-pair density and the current density distribution. 
By putting a smaller superconducting disk or ring in the center of a larger ring, the properties change 
drastically. Extra ground state transitions are found, where the total vorticity stays the same, but 
the vorticity of the inner superconductor changes by one. Due to the magnetic coupling, the current 
in the external ring exhibits extra jumps at the transition fields where the vorticity of the inner 
superconductor changes. In this case, for certain temperatures, re-entrant behavior and switching 
on and off of the superconducting behavior of the rings are found as function of the magnetic field. 
A H — T phase diagram is obtained for the situation where the inner ring has a higher critical 
temperature than the outer ring. An analytic expression for the magnetic coupling is obtained for 
thin rings and extreme type-II superconductors. 

PACS numbers: 74.20.De, 74.25.Ha, 74.25.Dw 



I. INTRODUCTION 



Single mesoscopic superconducting ringsEFO and diskalj~E3 have attracted a lot of attention during the last years. 
The superconducting properties of these single samples have been studied experimentally and are rather well described 
and understood theoretically. ___ _ __ 

Many of these studies were limited to disksEjl3E!J~Ej and ringsD'Ll of zero thickness. In such a case the magnetic field 
induced by the supercurrents are neglected and one assumes that the total field equals the external field, which is taken 
uniform. In the present paper disks and rings with finite thickness are considered. In this case the superconductor 
tries to expell the field. Supercurrents are induced creating a local magnetic field opposite to the external one. This 
leads to a decrease of magnetic field inside the superconductor and an increase near the sample edges. For sufficiently 
large disks and magnetic fields, the magnetic field penetrates the superconductor and vortices are created, which 
corresponds to regions of local field compression. In small disks, only circular symmetric vortex states nucleate in 
the center, which are called giant vortex states. Encircling the vortex on a closed loop, the phase change of the order 
parameter equals L times 2ir, where L is the angular momentum or vorticity. In larger disks, seperate vortices can 
nucleate creating a multivortex state and the vorticity L is just the number of vortices. In the case of rings, the 
magnetic field can be expelled to the inside or to the outside of the ring or partially to both sides. 

In the present paper we want to understand what will happen if two mesoscopic superconductors are put close to 
each other but are electrically isolated. When an uniform magnetic field is applied, it will be locally altered by each 
of the superconductors. In some regions of space the field will be expelled from the superconducting disk or ring, 
while in other regions it will be compressed into vortices penetrating the sample or compressed into the inside of a 
superconducting ring. This results in a strongly nonuniform total field, which is the superposition of the applied field 
and the field created by the supercurrents. Another superconductor will interact with this nonuniform field. In this 
way the superconductors are coupled by the magnetic field and they are interacting with each other. This coupling 
will influence the properties of the two superconductors. It is this coupling that will be studied in the present paper. 

The magnetic coupling between normal mesoscopic rings is well known. Wang et ac3 studied theoretically the 
persistent currents of N normal rings placed periodically on the same plane. Because of the mutual inductance 
between two rings, the electric current in one ring produces an induced flux in the other ring, creating an extra 
current in this ring. They found that the mutual inductance between normal rings enhances the persistent current as 
was the case in the experiments described in Ref£3. 

Correlations in strays of magnetically coupled superconducting aluminium rings were investigated experimentally 
by Davidovic et a£a They used ultrasensitive susceptibility techniques and scanning Hall probe microscopy to study 
arrays of electrically isolated superconducting rings of micron-size. When the external flux is close to 3>o/2, the 
magnetic moments produced by the supercurrents in such rings are analogous to Ising spins. Magnetic moments 
parallel to the applied magnetic field can be called "spin up", while those in the opposite direction "spin down". 
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Via their dipolar magnetic fields, neighbouring rings can interact antiferromagnetically and the different rings are 
influencing each other. ■— ■ 

Recently, Morelle et a£Q studied experimentally the magnetic interaction between two superconducting mesoscopic 
aluminium rings, close to the superconducting/normal phase transition. In their sample, a smaller ring was placed 
in the center of the larger ring. Using resistivity measurements the phase boundary was obtained for the two-ring 
structure as well as for the reference single ring. In both systems, Little-Parks oscillations were observed in the H — T 
phase diagram. The modification of the T C (H) oscillations of the outer ring is seen in the Fourier spectrum of the 
T C (H) line due to the coupling between the outer and the inner ring. They suggested that an inner ring made from 
a different superconductor with a higher critical temperature would increase the magnetic coupling between the two 
rings. 

In the present paper we present a theoretical investigation of the properties of two coupled mesoscopic supercon- 
ductors. Our main attention will go to the interaction between the two superconductors. How are they influencing 
each other? How do the superconducting properties of a single ring change when another superconductor is placed 
in its center? Therefore, we consider two different configurations: (i) a ring-disk configuration where a small disk is 
placed in the center of a larger ring, and-Zii) a ring-ring configuration where a small ring is placed in the center of the 
larger ring as in the experiment of Ref.c3. We will also give an example of a ring-ring system where the inner ring is 
made from a different superconductor with a higher critical temperature. Our theoretical analysis is based on a full 
self-consistent numerical solution of the coupled non-linear Ginzburg-Landau equations. Since we consider sufficiently 
narrow rings and small disks, only axial symmetric giant vortex states will nucleated. Therefore, the equations can be 
solved for a fixed value of the vorticity. The magnetic field profile outside and inside the superconductor is obtained 
self-consistently and, therefore, the full demagnetization effect is included in our approach. 

The paper is organized as follows: In Sec. II we describe the theoretical formalism. Our results for the ring-disk 
configuration are discussed in Sec. III. In Sec. IV the ring-ring configuration is studied. In Sec. V we calculate the 
H — T phase diagram for the two ring system, where the inner ring is made of a different material with a higher 
critical temperature. In Sec. VI we analytically calculate the energy of two coupled thin rings made of a different 
material for high values of k. Finally, in Sec. VII our results are summarized. 



II. THEORETICAL FORMALISM 



We consider a superconducting ring with inner radius Ri , outer radius R and thickness d immersed in a insulating 
medium (for example vacuum). In the center of this ring a superconducting disk with radius R* or another super- 
conducting ring with inner radius R* and outer radius R* is placed with the same thickness (see Fig. [I]). The whole 
sample is placed in a perpendicular uniform magnetic field M — (0,0, Ho). To solve this problem, we expand our 
previous approach for thin superconducting disksLJ to a system of two axial symmetric superconductors each made of 
a different material. In the present paper we solve the system of two coupled nonlinear Ginzburg-Landau equations 
which determine the distribution of both the superconducting order parameter \P ( r ) and the vector potential A ( r ) 
inside and outside the superconductor, 



ttf- — * = -a il0 *-A, o *|*| 2 , (la) 




VxVx^ = —y, (ib) 

c 

where for the inner superconductor the parameters are ct; j0 = ct;, /3j j0 = /?;, and for the outer ring a,^ — a o, Pi,o — O - 
The density of superconducting current j is given by 

J* = — ( ^1- - W^*) - — |*| 2 1 , (2) 
im V / mc 

and m is the mass of the Cooper-pair. Since we only consider circular symmetric rings and disks, we use cylindrical 
coordinates. Any position in space will be expressed by ~r* = (p, </>, z) , where p is the radial distance from the center, 
<f> is the azimuthal angle and z is the perpendicular direction. The sample lies between z = —d/2 and z = d/2. 

Eqs. (|la|)-(||) have to be supplemented by boundary conditions for ^ (T*) and A (T*). The boundary condition for 
the superconducting condensate at the superconductor/insulator surfaces is given by 
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(3) 



p=R* 1 ,R*,R 1 ,R 2 



where ~rt is the unit vector in the radial direction. This condition expresses that the superconducting current in the 
radial direction vanishes at the sample surface. The boundary for the vector potential has to be taken far away from 
the disk 
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where the field equals the applied magnetic field Tt = (0, 0, Hq). 

To rewrite Eqs. (|lal)-(p|) in dimensionless variables, we express all distances in units of the coherence length of the 
outer superconductor £ D = %/ y/~2ma , the order parameter in ^o.o = y/— ot a / /3 , the vector potential in ch/2e(; 
and the magnetic field in H c i,o = c ^/2e£ 2 = n y/2H C0 , where H co = yj Ana 2 / ' (3 is the thermodynamical critical 
field of the outer superconductor and k = A /£ is the Ginzburg-Landau parameter for this superconductor. The 
penetration depth of the outer ring is given by A D = 
London gauge, div A = 0, we can rewrite Eqs. (lis 



y/m/n (c/4e) /^o,o- Using these dimensionless variables and the 
-(0) in the following form: 
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where £j jG = £j, Ki j0 = Ki in the inner superconductor and £j )0 = £„, Ki >0 = k in the outer ring. The first boundary 
condition becomes 
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p=R{,R%,Ri,R 2 



* = 
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Provided that Eqs. jla|)-(p|) are fulfilled, the difference between the Gibbs free energy of the superconducting state 
and the normal state is determined by the expression 



'f m' 



(7) 



where the integral is over the total volume V of the superconducting samples and A$ = ^it^Hop is the external 
vector potential in the absence of a superconductor. The free energy is measured in H 2 V/8ir. The dimensionless 

supercurrent ~~j is given by ~j = HP*^* - * V***) /2i - |*| 2 ~A. 

We restrict ourselves to samples with thickness d < £, which implies that we are allowed to assume that the order 
parameter does not vary in the z directio rOO On the other hand for the vectoc-potential the variation in the z 
direction is retained and for d > A this is very important due to the Meissner effectcjO. 

We consider sufficiently narrow rings and small disks, and therefore only axial symmetric giant vortex states will 
nucleateoO. Consequently, the equations can be solved for a fixed value of the angular momentum L out in the outer 
ring and Li n in the inner superconductor that leads to the order parameter 

where Li n ,out = L ou t in the outer ring and Lin, out = Li n in the inner superconductor. Consequently, both the vector 
potential and the superconducting current density are directed along the azimuthal direction it^. For fixed angular 
momenta L out and Li n , Eqs. (|5a|)- (|5b|) can be reduced to the following form: 



-ld_ 
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(8b) 
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where 9(x) = 1 for x < 1 and 8(x) = for x > 1, A = A~$ $ and () indicates averaging over the disk thickness 

Because the superconducting condensates of the inner and outer superconductors are disconnected from each other 
they can not influence each other directly. The coupling is entirely due to the magnetic field, or equivalently the 
vector potential, as is expressed by the second term in Eq. (|8a|). The total magnetic field is a sum of the applied 
field and the field created by the superconducting currents of the inner ring and the outer superconductor which is 
described by Eq. where / (p) = fi„ (p) + f ou t (p) and /,„ (p) [f ut (p)\ is only different from zero in the interval 
R* <p< R* [Ri < p < R ] . 

The magnetic field created by the supercurrent has a H ~ 1/r 3 dependence for large r. Therefore, we can restrict 
our calculations to a region with radial size R s which we took typically 5 times the sample size and longitudinal size 
d s typically 10 times the sample thickness . The boundary conditions for the outer parameter can be written as 



df 
dp 



= 



(9) 



p=R*,R*,R 1 ,R 2 



for the ring-ring configuration, and 



dp 



= > P 7T 

p=R',R 1 ,R 2 °P 



, 



(10) 



p=0 



for the ring-disk configuration. The condition 
boundaries of our simulation region 



for the vector potential taken at infinity can be transferred to the 



A{z,p: 



Rs) — ^H R S 



A(\z\ = d s ,p) = ^H Q p 



(11) 



Following the approach of Ref.0 we apply a finite difference representation on the space grid p n ,z m to solve 
Eqs. (pa|)-(Sb). Since the size of our simulation region exceeds by far those of the sample, we apply non- uniform 



space grids to diminish the computer time. The space grid is taken uniform inside the sample, and we increase the 
grid spacing exponentially with distance outside the sample. This allows us to use the same number of grid points, 
typically 128, inside and outside the sample. To obtain steady-state solutions for a system of two superconductors, 
the following iteration procedure was used: 
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\ Pn / 



where A m ^ n = A(z m ,p n ), f n = f (p n ), Pn+1/2 = (pn+1 +p n )/2 and z m+1/2 = (z m +i + z m ) /2. The upper index k 
denotes the iteration step. To speed up the convergency we introduced the iteration parameters 77/ and T] a , and we 

expanded the nonlinear term (/^) = (ft 1 ) 3 + 3 (ft 1 ) 2 (fn ~ ft 1 ) ■ 



III. RING-DISK CONFIGURATION 



First, we consider a superconducting ring with a superconducting disk in the center. We investigate the influence 
of the disk on the properties of the ring. As an example, we take a ring with inner radius Ri = 1.5£ and outer 
radius R = 2.0£ and a disk in the center with radius R* — 1.0£. Both superconductors have the same thickness, 
d = 0.15£, and thepGinzburg-Landau parameter was taken to be k = 0.28, which is a typical value for mesoscopic Al 
superconductorst-MLa. 
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Fig. g(a) shows the free energy for the considered system as a function of the applied magnetic field. First we 
consider the uncoupled system and calculated the free energy for the disk in the center (thick dashed curve) and the 
free energy for the different giant vortex states in thexuiter ring (thick dotted curve). The results for single rings and 
disks were exhaustively described in previous paperstZlu. Notice that, in this paper, the free energy is expressed in 
units of Fq — H^V/8tt, where V is the sum of the disk and the ring volume. This is-jLhe reason why the free energy of 
the disk and the ring are not equal to — Fq at zero magnetic field as it was in Refs.EHtl The size of the disk is such that 
only the Meissner state, i.e. the L,„ = state, can nucleate. At applied magnetic fields Hq/ H C 2 > 2.9 the disk is in the 
normal state, which results in F = 0. In the single ring, on the other hand, different giant vortex states with vorticity 
L ut = up to L out — 10 can nucleate before the ring becomes normal at Hq/H C 2 ~ 6.8. Next, we introduced the 
magnetic coupling between the disk and the ring and the results are given by the solid curves in Fig. |(a). The different 
axial symmetric states are determined by the vorticity of the disk L in and the total vorticity L out , which is equal to 
the vorticity of the ring. Therefore, we characterize the states by {L outl Li n ). For the considered configuration, we 
find states with Li n = and L out — up to L out = 5. We also find states with Li n = 1 and L out = up to L out = 10, 
which equal the giant vortex states of the single ring, because the disk is now in the normal state. Notice further that 
we could also write (6, 0) instead of (6, 1) because for the applied magnetic fields where the L ou t = 6 state in the ring 
exists the disk is normal even for L{ n — 0. We have chosen to write h\ n = 1 because this expresses more clearly that 
there is flux going through the disk. If both the disk and the ring are superconducting, the free energy of the total 
sample is different from the sum of the free energies of the single disk and the single ring. 

To investigate these new states in more detail we consider as an example the (2, 0) state. Figs. ||(a-c) show the 
magnetic field distribution, the current density and the Cooper-pair density, respectively, as a function of the radial 
position for five different applied magnetic fields, i.e., H /H c2 = 0.1, 0.5, 1.5, 2.0, and 2.5. Near H /H c2 = the 
(2, 0) state equals the L in = state of the disk and the ring is in the normal state. The reason is that the applied 
field is so low that a lot of magnetic flux has to be attracted to create a state with L out = 2 in the outer ring. 
Therefore, a very high superconducting current has to flow through the outer ring which leads to the destruction of 
superconductivity in this ring. The solid curves in Figs. ||(a-c) show that the Cooper-pair density and the current 
density are indeed zero in the ring. The magnetic field distribution shows the flux expulsion from the disk. Inside the 
disk the field decreases and at the edge there is a peak which illustrates a higher concentration of field because of the 
demagnetization effects. With increasing external field less flux has to be attracted and the current in the outer ring 
decreases. At Hq/H C 2 ~ 0.17 superconductivity is restored in the external ring (see the dotted curves in Fig. ^(c)). 
The dotted curves in Fig. ||(b) show that the current in the outer ring flows in the opposite direction than the current 
in the disk. The superconducting currents in the disk expel the flux, while the currents in the ring are attracting flux, 
which is compressed in the region between the disk and the ring (see the dotted curve in Fig. ||(a)). The free energy 
becomes now more negative as compared to the free energy of the single disk (see Fig. ^(a)). Increasing the magnetic 
field further leads to less attraction of flux and, hence, to a higher Cooper-pair density in the ring and a more negative 
free energy. When the external flux becomes comparable with the flux needed for the L out = 2 state, the outer part of 
the ring expells the flux to the outside, while the inner part of the ring still expells flux to the hole region. Therefore, 
the superconducting current changes sign in the ring region (see the dashed curve in Fig. |(b)). Since the flux is 
expelled in both directions, the dashed curve in Fig. ||(a) shows a positive peak at both ring boundaries. Further 
increasing the external field leads to external fluxes larger than the flux needed for the L out = 2 state and, hence, the 
ring has to expel flux in order to keep vorticity L out = 2. As a consequence, the current in the ring has to flow in 
the same direction as the current in the disk, which is shown by the dash-dotted curve in Fig. ||(b). Because of the 
expulsion, the field between the two superconductors is lower than the external field (see the dash-dotted curve in 
Fig. §(a)). If we further increase the magnetic field, the superconducting currents in the outer ring have to increase in 
order to expel more flux and consequently the Cooper-pair density decreases in the outer ring. At Hq/H C 2 ~ 2.4 the 
supercurrent becomes too high and the ring becomes normal again (see the dash-dot-dotted curves in Figs. ||(a-c)). 
At this field, the free energy equals the free energy of the single disk. 

The above discussion shows clearly the interplay between the superconducting state of the disk and the ring. Next, 
we investigate the interaction between the disk and the ring. Therefore, we added the sum of the free energies of 
the single disk and the single ring (thin dashed curves in Fig. ||(a)) and compare this result with the result of the 
ring-disk configuration (solid curves in Fig. ||(a)). Notice that there is a small difference between the two set of curves, 
which is due to the coupling between the two superconductors. Below we will show that this difference becomes more 
pronounced for thicker samples. 

Now, we will determine the attraction or expulsion of the magnetic field by the coupled superconducting system. 
Fig. ||(b) shows the magnetic field expelled from the sample, —M, as a function of the applied magnetic field: 

{H)--Ho 

where (if) is the magnetic field averaged over the area p < R Q , i.e. the outer radius of the ring and Hq is the applied 
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field. The thick dashed and thick dotted curves are the results for the single disk and the single ring and the solid 
curves for the total ring-disk system. By putting a disk in the center of the ring, more field is expelled and less field 
is attracted. Of course, for Hq/H C 2 > 2.9 the disk is in the normal state and we recover the result for the single ring 
case. 



IV. RING-RING CONFIGURATION 

In this section we replace the disk in the center by a second ring and the influence of this inner ring on the outer 
ring will be investigated. As an example for this ring-ring configuration, we consider a superconducting ring with 
inner radius Ri — 1.5£ and outer radius R = 2.0£ and a second ring in the center with inner radius R* = 0.6£ and 
outer radius R* = 1.1£. Both rings have the same thickness, d = 0.15£ and the same Ginzburg-Landau parameter 
k = 0.28. 

Fig. |] shows the free energy as a function of the applied magnetic field for the small single ring by thick dashed 
curves, for the larger single ring by thick dotted curves and for the coupled ring-ring situation by the thin solid curves, 
where the interaction between the two rings is taken into account. In the inner ring superconducting states can nucleate 
with vorticity Lj„ = 0, 1 and 2 and the superconducting/normal transition field is at Hq/H C 2 ~ 6.4. In the outer 
ring states with vorticity L out = up to L out = 10 exist and superconductivity is destroyed at H$/H C 2 ~ 6.75. The 
superconducting states nucleating in the double ring system can be characterized again by the indices (L out , L in ). For 
Lin = 0, superconducting states can nucleate with L out = up to 4, for Lj n = 1 with L out — 1 up to L out = 8, and 
for Li n = 2 with L out = 5 up to 10. For Li n 3 the states equal the states of the single outer ring because the inner 
ring will be normal. 

The indices in the figure correspond to the ground state of the coupled ring system. For the numerical example 
studied in the previous section the number and the position of the ground state transitions are the same as for the 
outer ring. In the present two ring system this is no longer the case and the number of transitions in the coupled 
system are larger than for the single outer ring case. The inner ring induces extra transitions in the coupled system 
each time when the vorticity of the inner ring changes with one unit. The first extra transition is the transition from 
(2, 0) to (2, 1) and the second one is the transition from (6, 1) to (6, 2). The (2, 1) state is the ground state in the 
magnetic field region 1.43 < Hq/H C 2 < 1.63 and the (6,2) state in the region 4.27 < H^/H C 2 < 4.28. Hence, by 
putting a ring in the center of the largeuring, the ground state shows extra transitions. This result corresponds to 
the experimental result of Morelle et aBB, who saw modifications of the T C (H) oscillations of the outer ring in the 
Fourier spectrum of the T C (H) line due to the coupling between the outer and the inner ring. 

Next, we focus further on the interaction between the inner ring and the outer ring. Therefore, we plot in Fig. |B| 
the ground state free energy of the coupled rings (solid curves) and the sum of the free energies of the two single 
rings (dashed curves) for the previous configuration (upper curves which are shifted by +0.1), i.e. d = 0.15£, and for 
a thicker sample with d — 1.0£ (lower curves). For d — 0.15£ the difference is most pronounced for the (2,0), the 
(5, 1) and the (6, 1) state. Both the value of the free energy and the transition magnetic fields are influenced by the 
interaction between the two rings. The left inset shows the (6, 1) — > (6, 2) — > (7, 2) transition in more detail. Notice 
that the interaction significantly decreases the magnetic field region where the (6, 2) state is the ground state. For 
d = 1.0£ the demagnetization effects become more important and therefore the interaction between the two rings gains 
importance. This results in a larger difference between the dashed and solid curves. The value of the free energy and 
the transition fields are changing considerably by the fact that both rings arc influencing each other. The two lower 
insets show the (2,0) — > (2,1) — * (3,1) and the (6,1) — ► (7,2) transitions in more detail. The magnetic field region 
where the ground state is given by the (2, 1) state decreases due to the interaction. For the (6, 1) — + (7, 2) transition 
coupling between the two rings leads to the interesting result that the (6, 2) state is no longer a ground state. 

Figs. ^|(a,b) shows the magnetic field expelled from the sample, —M, as a function of the applied magnetic field 
for the single outer ring (dashed curve) and for the double ring configuration (solid curve) for thickness d = 0.15£ 
and d = 1.0£ respectively. For d = 0.15£, the two extra transitions, resulting from the influence of the inner ring, are 
clearly visible by the jumps at Hq/H C 2 — 1.43 and 4.28. Notice that, depending on the direction of the current a single 
ring expels or attracts field at a given applied magnetic field and vorticity. Expulsion (attraction) leads to a lower 
(higher) magnetic field density in the center of the ring and to a higher (lower) density near the outside. Therefore, 
depending on the applied field, the field expulsion or attraction in the coupled ring configuration can either increase 
if the two currents are in the same direction, or decrease if the directions are opposite. This can be clearly seen from 
Figs. ^(a,b) where the field expulsion or attraction becomes more pronounced by putting the inner ring in the center 
of the outer ring. For d = 1.0£ one extra transition results from the influence of the inner ring, i.e. at Hq/H C 2 = 1.58. 
The difference between the transition fields of the single ring and the coupled ring system becomes larger (see also 
the lower curve in Fig. f*) and the difference in the expulsion becomes more pronounced with increasing the sample 
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thickness, which indicates again that the interaction between the two rings increases with increasing d. 

In Figs. @(a,b) the magnetic field range, A7?o, over which the (L out , Li n ) state is the ground state, is plotted as 
a function of L out for thickness d — 0.15£ and d = 1.0£ respectively. This magnetic field range corresponds to the 
distance between two consecutive jumps in the expelled field (see Fig. ^). The results for the single outer ring are 
given by the open squares and for the double ring system by the closed circles. The curves are guides to the eye. 
For d = 0.15£ the extra transitions are clearly visible at L out = 2 and 6 and for d — 1.0£ at L out — 2 where for 
the same L ou t two jumps occur due to a transition of the inner ring. Also for the other vorticities L ou t, there is a 
difference between AHq for the single ring and the double ring. The reason is that the ground state transition fields 
are influenced by the interaction between the two rings. This was also visible in Fig. ||. If the free energy of the double 
ring was just the sum of the free energies of the two single rings, AHq would be the same for the single outer ring and 
the double ring, except for L out — 2 and 6, where extra transitions occur because L in changes with one unit. Notice 
further that the difference between the results for the single outer ring and the double ring enhances with increasing 
sample thickness. 

Next, we investigate the effect of the interaction between the two rings on the superconducting current density in 
the two rings for d — 0.15£. Figs. ||(a,b) show the averaged current density for the ground state in the inner ring and 
the outer ring, respectively, as a function of the applied magnetic field. The results for the single ring are given by 
dashed curves, these for the double ring configuration by solid curves. First, we describe what happens if there is no 
interaction between the two rings. In this case we can consider them as two single rings. At low fields, the ground 
state of a single ring is given by the L = state or Meissner state and the ring expels the field to the outside of the 
sample. With increasing external field, more flux has to be expelled from the ring which leads to a higher current 
density. After the first transition the ground state is given by the L = 1 state and initially flux will be trapped in the 
ring and the flux going through the ring is larger than the flux of the external field. To compress this extra magnetic 
field, the superconducting current in the ring has to flow in the opposite direction. At the transition, the current 
shows a jump from a negative to a positive value, i.e. from expulsion to compression. With increasing external field, 
less flux has to be compressed to achieve vorticity one and the current density in the outer ring decreases. Further 
increasing the field, the external flux becomes larger than the flux needed for L — 1 and flux has to be expelled. 
Therefore, the current in the ring changes sign. Without interaction between the two rings, the current density in one 
ring exhibits only jumps when the vorticity of the ground state of this ring changes (see dashed curves in Figs. ||(a,b)). 

In the coupled two rings situation (j) shows small jumps on top of the previously described expulsion — > compression 
jumps. At low fields, the ground state is given by the (0, 0) state or Meissner state. Both rings expel the field to 
the outside of the sample, which means that the current flows in the same direction in each ring. Since some flux is 
already expelled by the outer ring, the inner ring has to expel less and therefore the current is less negative. After 
the first transition the groundstate changes into the (1, 0) state. Now, the outer ring compresses the field to achieve 
vorticity one, and, as a consequence, the field in the hole of the outer ring is larger than the external field. This means 
that the inner ring has to expel more field and the current density jumps to a value more negative than its value 
without interaction. The other transitions can be explained analogously. From Figs. ||(a,b) it is clear that the two 
rings are influencing each other and that the interaction between the two rings results in extra jumps in the current 
density in one ring when the vorticity of the other ring changes. These jumps are smaller than the jumps when the 
vorticity of the considered ring increases, but they are not negligible. 

Up to now, we considered rather small samples. For the single ring it is known that by increasing the sample size 
(i) more L states are possible and (ii) the magnetic field range over which the state with vorticity L is the ground 
state decreases! Therefore, for a larger radius of the double ring configuration, we expect many more ground state 
transitions. Fig. || shows the ground state free energy for a single inner ring with radii R* = 2.0£ and R* = 1.5£, 
for a single outer ring with radii R a = 3.0£ and Ri — 2.6£ and for the coupled ring-ring configuration. The sample 
thickness is d = 0.15£ and the Ginzburg-Landau parameter k — 0.28. For the single inner ring, the ground state 
changes from vorticity Lj„ = up to L{ n — 10 and the superconducting/normal transition field is at Hq/H C 2 = 6.73. 
For the single outer ring the ground state changes from vorticity L out — up to L out — 32 and superconductivity 
is destroyed at Ho/H C 2 = 8.40. By comparing the free energy of the double ring with the one of the outer ring, we 
notice that there are many more ground state transitions as a consequence of the transitions in the inner ring. For 
the single ring the minimum in the free energy of the L + 1 state is always less negative than the one of the L state. 
Due to the interaction between the two rings, this is no longer always the case for the double ring configuration. At 
Hq/H C 2 > 6.73 the free energy of the double ring configuration equals the one of the outer ring since the inner ring 
is in the normal state. 

Figs. |l(](a,b) show the field expelled from the region p < R* of the inner ring with radii R* = 2.0£ and R* = 1.5£ 
and the region p < R a of the outer ring with radii R a = 3.0£ and Ri = 2.6£, respectively. The results for single rings 
are given by the dotted curves and for the double ring configuration by the solid curves. At low fields, the single 
inner ring is in the Meissner state and expels the magnetic field, i.e. — M > 0. With increasing external field, more 
field is expelled. At H /H c2 ra 0.33 the ground state changes from L in = to L in — 1 and flux has to be compressed 
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into the hole to achieve vorticity one. Therefore, the magnetic field inside the hole will be larger than the external 
one and — M jumps to negative values, i.e. field compression. With increasing field, less flux has to be attracted and 
— M becomes less negative. Further increasing the external flux becomes larger than the one needed for Lj„ = 1 and, 
therefore, the field has to be expelled again. This means that —M becomes positive. Further increasing the field, 
more flux has to be expelled and —M becomes more positive. At Ho/H C 2 ~ 0.99 the vorticity changes from Lj„ = 1 
to Li n — 2, which means that — M jumps to negative values, and so forth. At H$/H C 2 = 6.73 superconductivity is 
destroyed and the field becomes equal to the external one, and as a consequence, —M = 0. The description for the 
single outer ring is completely analogous. 

Placing a larger ring around the inner ring influences the expelled field drastically (see the solid curves in Fig. |l0|(a)). 
Due to the expulsion of the outer ring at low fields, the magnetic field inside the hole of this ring will be smaller than 
the external one. Now, the expulsion by the inner ring results in a smaller local field than for the case of the single 
ring and, as a consequence, the expelled field increases at low fields. At H^/H C 2 ~ 0.13 the ground state of the outer 
ring changes from vorticity L out = to L out = 1, which means that suddenly the outer ring has to attract flux to 
achieve vorticity one. Therefore, the field inside the hole of the outer superconductor becomes larger than the external 
one and the expulsion by the inner ring will be less pronounced than for the case of the single ring. As a consequence, 
— M jumps from a value above the one for the single ring case to a value below this value. This interplay between the 
two rings leads to a higher expulsion (attraction) from the region p < R* when the outer ring expels (attracts) flux 
and to a lower expulsion when the outer ring attracts (expels) flux. For the outer ring an analogous explanation can 



Until now, we considered always two superconductors made of the same material. This means that both super- 
conductors have the same coherence length, penetration depth and critical temperature, i.e. = £ Q , Xi — A G , and 
T c ,i — T c , . Since both rings have the same width and the radius of the inner ring is smaller than the one of the 
outer ring, the inner ring becomes normal at a smaller field than the outer ring. As a consequence, no effect of the 
magnetic coupling can be observed in the H — T phase diagram. To circumvent this problem the T c of the outer 
ring was artificially lowered in the experiment of RcfE3 by applying a sufficiently large external current through the 
outer ring. An alternative approach will be followed in the present section where we take the inner ring of a different 
material such that it has a higher critical temperature than the outer ring, and also a different coherence length and 
penetration depth, which leads to a different Ginzburg-Landau parameter. _ . 

As an example, we take for the outer ring the values used by Geim et at3 for Al, i.e. £ Q (T = 0) = 250nm, 
A (T = 0) = 70nm and thus n = 0.28, resulting in a critical temperature T Ct0 (H = 0) = 1.3!^. For the inner ring, 
we assume a higher critical temperature T c ^ — l.2T Cj0 = 1.56-fT, and (T = 0) = 160nm, Xi (T = 0) = 80nm and 
thus Ki = 0.5. For the radii of the outer ring we take as an example Ri — 375nm, R a — 500nm, and for the inner 
ring R* = \2hnm and R* — 250nm. The H — T phase diagram is shown in Fig. for the uncoupled situation 
for the inner ring (thick dashed curves) and for the outer ring (thick dotted curves) and the coupled double ring 
situation (solid curve). At T = the outer ring has a much higher critical field (H nuc / H C 2. — 6.74) than the inner 
ring (H nuc / H C 2_ = 4.34). Therefore, the superconducting/normal transition of the double ring configuration equals 
the one of the outer ring for low temperatures. With increasing temperature, the nucleation field of the outer ring, i.e. 
the one of the double ring system, changes more quickly than the one of the inner ring. The oscillations are the well 
known Little-Parks oscillations. At T/T c _ a — 0.912 both single rings have the same transition field H nuc / H C 2 i0 = 1.91. 
At higher temperatures, the superconducting/normal transition is determined by the inner ring. 

The situation where the critical field of the outer ring is larger than the one of the inner ring is exhaustively 
described in the previous sections. In Fig. [l^ we show the free energy for the configuration of Fig. O at T = 0.98T CjO 
where the superconductivity of the inner ring exists at larger fields than the one of the outer ring. The free energy of 
the superconducting states of the inner ring are given by the dashed curves, the states of the outer ring by the dotted 
curves and the double ring configuration by the solid curves. Both in the inner and the outer ring, superconducting 
states with vorticity L — and L — 1 exist. At T = 0.98T C O , the critical fields of the inner and the outer ring are 
Hq/H C 2, = 1-77 and 0.69, respectively. Notice that in both rings the free energies of the L = state and the L = 1 
state do not cross, which means that with increasing field the ground state changes from the Meissner state into the 
normal state and, with further increasing the field, into the L = 1 state and back into the normal state. The reason 
is that near T c the superconductivity of the ring has decreased. This means that only rather small currents can be 
induced and thus only a small flux can be attracted or expelled by the ring. In the region between the existence of the 
L = state and the L = 1 state the currents, which have to be induced to expel or attract the necessary flux to achieve 
vorticity L — and L = 1, are too high. With increasing temperature, the L = 1 state can not nucleate anymore 
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and the superconducting/normal transition jumps to the field where the L = state is destroyed. The corresponding 
oscillations in the H — T phase diagram (Fig. [by) are the Little-Parks oscillations. For the double ring configuration 
the (0,0), the (1,0) and the (1, 1) state can nucleate. The (1, 1) state is split into two parts corresponding to the 
L = 1 states in the two single rings with an intermediate magnetic field region in which both superconductors are 
normal. The ground state changes from the Meissner state (0,0) into the (1,0) state at Hq/H C 2^ ~ 0.53, which 
equals the (1, 1) state at H /H C 2 vO > 0.69. Further increasing the field the ground state changes into the normal state 
at Hq/H C 2. ~ 0.76, then back into the (1,1) state at Hq/H C 2 i0 ~ 0.86 and further back into the normal state at 
Hq/H C 2, ~ 1.77. Compared to the uncoupled inner ring and the outer ring situation, extra ground state transitions 
occur for the double ring case with interesting re-entrant superconducting behavior and a switching on and off of the 
superconducting state in the inner and outer ring. 



VI. TWO COUPLED THIN RINGS IN THE LIMIT k > 1 



Here we will show that in the limit of two coupled thin rings it is possible to obtain analytical results for the 
coupling energy between the two rings. This also corresponds to the case of k ^> 1 and allows to solve the problem 
analytically with the small parameter d/ n 2 <C 1. From the numerical calculations of previous sections it follows that 
the radial dependence of the order parameter in both inner and outer rings is slow and smooth. Therefore, we assume 
that the order parameter in both inner and outer rings, fi n and f ou t, respectively, are constant. In the limit of the 
thin rings we neglect in the first approximation the z— dependence of the vector potential. With this assumption and 
because of the cylindrical symmetry of the problem the vector potential has only the azimuthal component A (p) and 
the magnetic field has only the normal component H (p) = d (pA (p)) /dp. 

The distribution of the vector potential due to the supercurreuts inside the inner and outer ring are described by 
the following equation (see also, for example, Eqs. (2)-(3) in Ref.EJ): 
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(13b) 



The corresponding boundary conditions are the continuity of both A (p) and H (p) on the radial sides of the rings. 
For d — > it gives everywhere A(d = 0) = H p/2 and H(d = 0) = H . To a first approximation the solution of 



Eq. (13a) becomes for small d/n 2 . 
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which is valid inside the superconductors, i.e. in the r ang e R* < p < R* and Ri < p < R . The vector potential 
outside the rings is then obtain from a solution of Eq. ( |13b| ) and is equal to 



A{ P )={ D 3P /2 + C 3 /p, 
H p/2 + C ext /p, 



o<p<R*, 
K<P< Ri, 
P>Ro- 



(14b) 



The corresponding distribution of the magnetic field is 

f 2 n {d/n 2 ){L m \np-Hop 2 /A) 
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(15) 



From the above boundary conditions the integration constants D and C in Eqs. (14a)-(|15|) are 







C' m = {d/An 2 )f 2 n R* 2 (-L ln + H Rf/A) , 
C 3 = (d/4 K f)ff n (Rf - Rf) [L in - H (Rf + Rf) /4] , 
C out = C 3 + (d/4n 2 )f 2 ut R 2 (-L out + H R 2 /A) , 
C ext = C 3 + (d/4K 2 a )f 2 at {R 2 - R 2 ) [L out - H (R 2 a + R 2 ) /4] , 
{d/K 2 )f 0Ut [L out lnR -HoR 2 /4], 
Ho + {d/K 2 )f 0Ut [L^ In (Jio/iZi) - HoS out /4iT} , 
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D 1= H + (d/K, 2 )f 2 n [L m In (R*JR*) - H Q S m /4n] + {d/ K 2 )f 2 ut [L out In (i? /i? 4 ) - H S out /4ir 



Dout 
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Notice that within this approximation the magnetic field is constant in the inner hole and in the space between the two 
rings, and has a small (as the ring's thickness is small) radial dependence inside the rings. Inserting the expressions 
for the order parameters, the magnetic field and the vector potential into the expression for the energy (note that in 
the equivalent expression (0) the Ginzburg-Landau equations have already been used) 
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With a O (d 2 /n 4 ) accuracy we find the difference between the Gibbs free energy of the superconducting state and the 
normal state 
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After the minimization of the free energy ( |l7| ) with respect to /i n ( ou t) we obtain two equations 
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from which we obtain the equilibrium values of the order parameters for the three possible situations. 

When both the inner and outer rings are in the superconducting state (case I), we obtain with a O (d 2 j 'k 4 ) accuracy 
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Inserting these expressions into Eq. ([l7]) we obtain F 1 = F* n + F^ ut 
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is the self energy of the inner ring, and 
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is the self energy of the outer ring, while 



2d ( k 2 




(21c) 



is the interaction energy between the two rings. 

When one of the rings is in the normal state the results become more simple. Namely, when the outer ring is in 
the normal state and the inner ring is superconducting (case II) or vice verca (case III) , we have 
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respectively. The corresponding energies are 
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and in the d/ -C 1 limit they coincide with i 7 ^ and -F/ ut , respectively. 

One can see, that an interaction between the two rings (i.e., the coupling) exists only when both rings are supercon- 
ducting. The energy of the ring-ring coupling in the considered limit is proportional to the ring's thickness. Due to 
the interaction between the rings the Cooper pair density in each ring has a small (proportional to d/ n 2 ) contribution 
from the neighbouring ring. 

In Fig. [l^ the calculated magnetic field dependence of the ground state energy of the coupled thin rings of the 
same material (i.e., with £j = £ D and ft, = k ) with the same radial sizes as in Fig. 5 is shown for d/n 2 = 0.05 (solid 
curve). Also for comparison the curves from Fig. 5 are shown by dashed and dash-dotted lines which correspond 
to d/n 2 ~ 1.98 and 12.76, respectively. One can see that all curves have the same qualitative behaviour and with 
increasing d/ n 2 only a small decrease of the ground state energy takes place. Also the values of the transition fields 
between the different L states are very nicely reproduced. This is quite surprising and, to make the physics of this 
result more clear, we consider the limit of narrow rings. 

Let's introduce the average radius of the two rings, p in — (R* + R*) /2,p out = (R + Ri) /2, respectively, and their 
corresponding widths, 2wt n — R* — R*,2w out — R — Ri. Next we expand Eqs. ( |l8| ) and (21) with respect to the 

small parameters 2wi n [ out )lJ> in ^ out ^ *C 1 with an O {w 2 n( ^ out ^j accuracy. Within this approximation the energies (21) 
are 
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where $i n (out) = HoPi n r out \/2 is the average magnetic flux through the corresponding ring, measured in units of 
$o = hc/2e. From these expressions one can see that the dominant thickness dependent terms in F?, ut y as well as 

the ring's interaction energy F* nt , arc of order O ((^/ K i( )) w in(out)^j ( n °t of order O (d/n,^ ^j, as one could think 
naively). This additional smallness shows that an increase of the ring thickness influences the ring's self-energies and 
the interaction energy very slightly and for rings which have a not too large width (as in our case) the above results 
become valid for d/n?^ > 1. 

VII. CONCLUSIONS 

We investigated the magnetic coupling between two concentric mesoscopic superconductors with non-zero thickness. 
When a second superconductor is placed in the center of a superconducting ring, it feels a non-uniform field, which 
is the superposition of the uniform applied field and the field expelled from the outer ring. Also the first ring will 
be influenced by the magnetic field expelled from the superconductor in the center. So, both superconductors are 
coupled magnetically. This results in substantial changes of the superconducting properties. 

From the study of the free energy we learned that extra ground state transitions occur in comparison with the single 
ring case. These are transitions where the total vorticity stays the same, but the vorticity of the inner superconductor 
changes by one unit. We also found that the free energy of the double ring system is not exactly the same as the 
sum of the free energies of the two uncoupled single rings which is another signature of the magnetic coupling of 
both rings. This interaction enhances with increasing sample thickness. We also calculated the expelled field for 
the ring-ring configuration which showed that as compared with a single ring more, or less, field can be expelled or 
attracted depending on the vorticities of both superconductors. 

The behaviour of the Cooper-pair density, the magnetic field profile and the current density was calculated. Since 
an extra superconductor is placed in the center, the magnetic field will be expelled from this superconductor or 
will be compressed in the center of it, which results in a higher or a lower magnetic field density between the two 
superconductors. The current in both rings exhibits extra jumps at the transition fields where the vorticity of the 
other ring increases or decreases by one. The reason is that at these applied fields the total magnetic field in the 
region between the two superconductors changes. 

We investigated what happens if the inner ring is made of a different material with a higher critical temperature. 
The H — T phase diagram showed that the nuclcation field of the double ring equals the one of the outer ring at low 
temperatures and the one of the inner ring at higher temperatures. 

Analytical expressions arc obtained for the magnetic field distribution and the energy of two coupled thin type-II 
superconducting rings. These analytical results are found to give excellent results when d/n?^ < 1 and, moreover, 
give also good agreement with our "exact" numerical results for sufficiently narrow rings when d/nf^ > 1. 
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FIG. 1. Schematical outline of the considered configurations; the ring-disk configuration (left) and the ring-ring configuration 
(right). 
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FIG. 2. (a) The free energy and (b) the expelled field as a function of the applied field for a ring with inner radius Ri — 1.5£ 
and outer radius R = 2.0£ (thick dotted curves) and a disk in the center with radius R* = 1.0£ (thick dashed curves) and for 
the coupled ring-dot configuration (thin solid curves). All superconductors have the same thickness, d = 0.15£, and k = 0.28. 
The thin dashed curves give the sum of the free energies of the single disk and the single ring. 
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FIG. 3. (a) The magnetic field distribution, (b) the current density and (c) the Cooper-pair density as a function of the 
radial position for the (L out , L in ) = (2, 0) state of the ring-dot configuration of Fig. | at H /H c2 = 0.1, 0.5, 1.5, 2.0, and 2.5. 
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FIG. 4. The free energy as a function of the applied magnetic field for a superconducting ring with radii Ri — 1.5£ and 
Ro = 2.0£ (dotted curves), a ring with R* = 0.6£ and R* = 1.1£ (dashed curves) and the double ring configuration (solid 
curves). Both rings have the same thickness, d = O.I5£, and n = 0.28. The indices indicate the ground state vorticities 
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FIG. 5. The ground state free energy of the double ring configuration of Fig. |4| (solid curves) and the sum of the free energies 
of the two rings (dashed curves) for d — 0.15£ (upper curves) and for d = l.Of (lower curves). The upper curve is shifted by 
+0.1 for clarity. The insets show some of the crossings in more detail. 
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FIG. 6. The magnetic field expelled from the sample, — M, as a function of the applied magnetic field for the single outer 
ring (dashed curve) and for the double ring configuration (solid curve). The sample is the same as in Fig. ^ with thickness 
d = 0.15£ (a) and d = 1.0£ (b). 
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FIG. 7. The magnetic field range, AHo, over which the (L ou t, Li n ) state is the ground state, as a function of L out for the 
single outer ring (the dashed curve and open squares) and for the double ring system (the solid curve and closed circles). The 
sample is the same as in Fig. H with thickness d = 0.15£ (a) and d — 1.0£ (b). 
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FIG. 8. The averaged current density for the ground state in the inner ring (a) and the outer ring (b) as a function of the 
applied magnetic field for the same double ring configuration as in Fig. ^. The results for the single rings are given by dashed 
curves, these for the double ring configuration by solid curves. 
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FIG. 9. The ground state free energy for a single inner ring with radii 7?* = 2.0£ and Rt = 1.5£, for a single outer ring 
with radii R = 3.0£ and R t = 2.6^ and for the double ring configuration, i.e. the combination of these two rings. The sample 
thickness is d = O.I5£ and the Ginzburg-Landau parameter k = 0.28. 
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FIG. 10. The field expelled from the region p < R„ of the inner ring with radii = 2.0£ and R* = 1.5£ (a), and the region 
p < R of the outer ring with radii R a = 3.0£ and Ri = 2.6£ (b). The results for single rings are given by the dotted curves 
and for the double ring configuration by the solid curves, (d — 0.15£ and k = 0.28). 
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FIG. 11. The H — T phase diagram for the inner ring (dashed curves), the outer ring (dotted curves) and the double ring 
configuration (solid curves). The material parameters and the sizes of both rings are different and are given in the figure. 
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FIG. 12. The free energy as a function of the applied magnetic field for the inner ring (dashed curves), the outer ring 
(dotted curves) and the double ring configuration (solid curves) for the system of Fig. |ll| at T — 0.98T c , o . The inset shows an 
enlargement of the phase diagram (Fig. |ll|) in the T/T c>0 ~ 1 region. 
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FIG. 13. The groud state free energy of the double ring configuration with the same parameters as in Fig. ^ obtained from 
our analytical expressions which are valid for type-II superconducting rings. Also the ground state energies from Fig. ^ are 
shown. 
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